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Abstract 

A series of Higgs field quantum numbers in the anti-grand unification 
model, based on the gauge group SMG"^ x ?7(1)/, is tested against the 
spectrum of quark and lepton masses and mixing angles. A more precise 
formulation of the statement that the couplings are assumed of order unity 
is given. It is found that the corrections coming from this more precise 
assumption do not contain factors of order of the number of colours, = 
3, as one could have feared. We also include a combinatorial correction 
factor, taking account of the distinct internal orderings within the chain 
Feynman diagrams in our statistical estimates. 

Strictly speaking our model predicts that the uncertainty in its pre- 
dictions and thus the accuracy of our fits should be ±60%. Many of the 
best fitting quantum numbers give a higher accuracy fit to the masses and 
mixing angles, although within the expected fiuctuations in a x^- This 
means that our fit is as good as it can possibly be. 



1 Introduction 



Almost the only window of information available to glimpse the laws of na- 
ture beyond the Standard Model is at present provided by the ca 20 parame- 
ters, mainly already measured but not predicted by the Standard Model itself. 
Among these, the 13 parameters comprising the masses, mixing angles and CP 
violating phase for the fermions — quarks and leptons — make up more than half. 
We have earlier put forward a model as a candidate for explaining the orders 
of magnitudes of the masses for the six quarks and the three charged leptons. 
In this anti-grand unification theory (AGUT) the smallncss of most of the Stan- 
dard Model Yukawa couplings is to be understood as a consequence of there 
existing some approximately conserved gauge charges forbidding these Yukawa 
couplings We actually managed in our AGUT model to successfully fit all 
the 9 charged particle masses]^ and the three mixing angles using only three 
Higgs field vacuum expectation values, in addition to the usual Weinberg Salam 
Higgs field. 

In the present article we want to put forward a few refinements of this model: 

1) First we want to take into account that the number of different Feynman 
diagrams that can contribute to a given mass matrix element — a given effective 
Higgs- Yukawa coupling — becomes rather large when the diagram is complicated, 
because it involves many Higgs field-lines. This number of diagrams can be 
expressed in terms of factorials and can easily be so large that, even for order 
of magnitude calculations, we should include them. 

2) In our previous work on this AGUT-fermion mass model, we introduced 
a Higgs field S for which we could take the vacuum expectation value (VEV) in 
the fundamental (Planck) units to be of order unity. Inclusion of the factorial 
corrections, just mentioned under point 1), actually makes the approximation 
of the S field VEV by unity give a poorer fit to the data. So, in order to regain 
as good a fit as before, we are forced to allow the VEV of the Higgs field S to be 
a variable parameter when the factorials are included; although its fitted value 
still turns out to be rather close to unity. 

3) If there are many types of particles in the model, or if we have particles 
that like quarks are say triplets under colour symmetry]^, one would a priori 
imagine that the number of particles or the number of colours — i.e. 3 — were 
tacitly treated as being of order unity in our previous fits of the quark and 
lepton masses. If this were really so, one would expect corrections of the order 
of, say, a factor 3 in our predictions. Actually our previous fits 0] always agreed 
with data to better than a factor 3, except for perhaps the most sensitive of our 
predictions — the CP-violation strength — which is especially uncertain, because 
it is a product of very many only order of magnitudewise known factors; in fact 
it was not included in the fit but simply predicted by our model. Even this 

^We do not discuss the neutrino mass problem here, which requires the introduction of 
another mass scale and corresponding Higgs particle. 

^In the AGUT model there is even a different colour group for each proto-generation. 
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deviation of our CP-violation prediction from the experimental number is only 
just about a factor 3. In order to make sense of such a good agreement for our 
previous fits, we have to formulate a more precise form of the loose statement 
that the coupling constants are all of order unity at the Planck scale; so that 
we work to an accuracy where the number of colours (for the three SU(3) gauge 
groups in the model) and the number of species of particles, with masses of order 
the Planck scale, become large compared to the "order of magnitude accuracy" 
to which we hope to fit. We must hope — and indeed the estimates in section 
6 suggest the hope is fulfilled — that we do not get large factors depending on 
the number of Planck scale species postulated in the model, because we do not 
know how many they should be. Rather we must show that their number is not 
important, provided they are at least there at all. 

It is the purpose of the present article to include the two first improvements 
1) and 2), and to discuss the third potential correction 3) which, however, turns 
out to actually only give somewhat milder corrections. For example we show 
that our prediction of order of magnitude degeneracy for quark and lepton 
masses in the same generation (except for the t and c quarks) remains valid, 
without any order Nc = 3 factors. So it becomes understandable that we can 
get order of magnitude fits with only about 40% deviations from the data. 

Once we embark on including factors of suppression connected with our 
Higgs field S, which were previously set equal to unity, it matters how many 
times the nonzero VEV of this field S is needed to generate a given mass matrix 
element (or effective Yukawa coupling). The required number of factors of this 
VEV in turn depends on the gauge quantum numbers of the other AGUT Higgs 
fields, which we called W , T and ^. So their quantum numbers must now be 
specified completely and not just modulo the quantum numbers of our Higgs 
field S, as was sufficient as long as the VEV of the S field could be considered to 
be of order unity. This specification of the precise gauge quantum numbers for 
the Higgs fields W, T, ^ and the field ^ws, containing the phenomenological 
Weinberg Salam Higgs field, gives rise to an infinite number of models. How- 
ever, one should first of all remark that it is extremely reasonable to assume 
that, in a specific version of the model, the charges can only take on discrete — 
essentially integer in appropriate units — values, and secondly we do not expect 
them to be very large. So, at the end, the number of models to be considered is 
rather limited. We shall label the various considered models by a set of integers 
{a, f3, 7, S) defined below in section ^. 

In the following section || we give a short review of the AGUT model. Then, 
in the next section |^, we explain and implement the factorial correction coming 
from point 1): the large numbers of permutations of the Higgs fields in the 
Feynman diagrams. In section ^ we then define the set of models obtained 
by the various specific quantum number choices for our Higgs fields, i.e. by 
choosing (a, /?, 7, <5). In section |^ we present the results of fitting quark and 
lepton masses and mixing angles with the various parameter choices and with 
the factorial corrections included. Also dicussed in this section are two variants 
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of how one interprets the order of unity factors in the calculation of the predicted 
masses etc.: one can either just put random phase factors on the products of the 
suppression factors and the factorial corrections, or one can further provide them 
with a random factor of order unity. The latter variation is not so important 
when one term dominates but when, as for the case of the d-quark mass in 
our model, more than one term contributes, it may give a correction. The best 
fitting Higgs field quantum number combinations are selected. Then follows, 
in section ||, our proposal for making the concept of the couplings being of 
order unity more precise. The outcome of this discussion is really that the 
results of section ^ are not significantly modified by the more careful definition 
of what the couplings being of order unity means. In section Q we present the 
argument that, since many Feynman diagrams contribute, each mass matrix 
element in our model should have a Gaussian distribution in the complex plane. 
This, in turn, has the significance that we can even calculate an expectation for 
what the random element in the masses should be and, thus, predict the degree 
of deviation between our fit and the experimental numbers. Indeed we shall 
see that our prediction for the uncertainty agrees well with the experimental 
deviation of the fit. Finally in section |^ we present our conclusions. 

2 Our AGUT-model, a short review 

Anti-grand unification is a well suited name for our AGUT model in a number 
of senses, as we will now explain. It is similar to the usual grand unified theories 
(GUTs), which are based on SU{5), S'O(IO) or some other simple group, in the 
sense that we postulate the extension of the Standard Model Group ||^ , SMG 
= S{U{2) X t/(3)) « f/(l) X SU{2) X SU(i), to a larger gauge group at a very 
high energy of order the Planck scale. However the AGUT gauge group can be 
considered "anti" , as it is not at all simple in the mathematical sense of the 
word. On the contrary, its Lie algebra consists altogether of ten cross product 
factors: four U{1) factors, three SU{2) factors and three SU{2>) factors. Below 
the AGUT scale there is no supersymmetry nor other new physics beyond the 
Standard Model. 

Further the usual GUTs unify by combining several Standard Model Group 
quark and lepton irreducible representations into the same irreducible represen- 
tation of the grand unified group. Our AGUT is "anti" in the sense that it is 
characterised by not uniting any of the Standard Model Group irreducible repre- 
sentations into a larger irreducible representation. In this way our model avoids 
the exact mass degeneracies predicted from the minimal SU{5) GUT, according 
to which the running masses of the quarks with charge —1/3 are degenerate 
with their charged lepton partners at the unification scale. Honestly speaking, 
except popssibly for the case of the b quark being degenerate with the r lepton, 
these predictions are wrong. In fact the unwanted predictions, rris = and 
nid = Tie, can only be tolerated by complicating the SU{5) model [Q. How- 
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ever it turns out that the AGUT model replaces these exact SU (5) predictions 
by only order of magnitude degeneracy predictions, which are compatible with 
experiment. 

While SU(5) is the smallest simple group unifying quarks and leptons, the 
AGUT group SMG^ x C/(l)/ is the biggest group, which keeps the Standard 
Model irreducible representations separate, in the following sense: there is no 
larger group, containing SMG^ x U{1) as a subgroup, which is faithfully repre- 
sented on the Standard Model quark and lepton fields alone (no right-handed 
neutrinos or other new fermions), without gauge or mixed anomalies. Indeed 
the SMG^ X U{\) group is uniquely specified as the group G extending the 
Standard Model group and satisfying the following four postulates [1). 

1. G should transform the three generations of Standard Model Weyl parti- 
cles into each other unitarily, so that G C [/(45). 

2. G should be anomaly free, even without using the Green-Schwarz 
anomaly cancellation mechanism. 

3. The fifteen irreducible representations of the Standard Model Weyl fields 
remain irreducible under G. 

4. G is the maximal group satisfying the other three postulates. 

Also the Standard Model gauge coupling constants do not unify in the AGUT 
model, but their values have been successfully predicted, at the Planck scale, 
using the so-called multiple point principle According to this principle, the 
coupling constants are fixed at such values as to ensure the existence of many 
vacuum states with the same energy density. We note that the top quark mass 
and the Weinberg Salam Higgs particle mass can also be predicted using this 
principle [||. 

Really the existence of this AGUT group means that, near the Planck scale, 
each of the three quark-lepton proto-generations has its own SMG gauge group 
and associated 12 gauge particles — i.e. the gauge bosons also come in genera- 
tions. In addition there is an extra abelian U(l) / gauge boson; the correspond- 
ing gauge charge Q/ is not carried by the left-handed fermions or any of the 
first proto-generation particles. So the non-zero U{1) / quantum numbers for 
the proto- quarks and leptons can be chosen as follows: 

Q/(TK) = Q/(&i?,) = Q/(ci?)-i (1) 
Q^{^^B) = Qf{sR) = Qf{tR) = -i. (2) 

We stress that the physical quarks and leptons are superpositions of the proto- 
fermions. In particular it turns out that the right-handed charm and top 
quarks are essentially permuted: the physical right-handed top quark is the 
Weyl component with second proto-gcneration quantum numbers, namely the 
right-handed proto-charm quark, to first approximation. 
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The representations of the Standard Model Group SMG = S{U{2) x C/(3)) 
satisfy the charge quantization rule: 

y/2 + d/2 + 1/3 = (mod 1) (3) 

Here y is the conventional weak hypercharge and d is the duality, which is 
defined to be when the weak isospin is integer and d = 1 when it is half 
integer. Similarly the triality t is defined to be for an SU{3) colour singlet 
and t — ±1 for a triplet or antitriplet respectively. In the AGUT model this 
charge quantisation rule is satisfied separately for the SMG quantum numbers 
associated with each quark-lepton generation. 

Now we turn to the breaking of the AGUT gauge symmetry group SMG^ x 
U{l)f down to the Standard Model group SMG, which is supposed to occur 
close to the Planck scale. Unlike for the quarks and leptons, the gauge quantum 
numbers of the Higgs fields responsible for this symmetry breaking are not de- 
termined by the theoretical structure of the model. We assume that these Higgs 
fields have "small" quantum numbers, belonging to singlet or fundamental rep- 
resentations of the non-Abelian groups. The charge quantisation rule, Eq. (||), 
then determines the non-Abelian quantum numbers from the Abelian ones. So 
it is sufficient to specify just the Abelian quantum numbers in the form of a 
U{1) charge vector: 

Q.(|,f.f,0,) ^ (4) 

where yt is the conventional weak hypercharge for the i'th proto-generation. 
For example the Higgs field introduced in |^] with a VEV of order unity in 
fundamental units, has the Abelian quantum numbers Qs — (i, — 0, — l)- It 
follows from the charge quantisation rule that the full set of AGUT quantum 
numbers for the S field are: (3, 2, i; 3, 2, — i; 1, 1, 0; — l) , where the sets of 
three quantum numbers specifying the representations under the three proto- 
generation S'AfG-groups are separated by semi-colons. 

As long as we consider the VEV of the Higgs field S to be of order unity, 
we are in reality working with the group to which S breaks the AGUT group 
down. Then it is only necessary to specify the quantum numbers under this 
subgroup (i.e. modulo the quantum numbers of S) . Thus the quantum numbers 
Qt — (Oj ~f ) for the Higgs field T used in could just as well be taken 

to be Qt — ^,^) modulo the S quantum numbers. This transformed 

set is of course arbitrary as far as adding S quantum numbers is concerned, 
but we have chosen it as in some sense the smallest Abelian quantum number 
choice. 

Another crucial assumption in our model is the existence, at the fundamental 
(Planck) scale, of a large spectrum of vector-like Dirac particles with the quan- 
tum numbers needed in the Feynman diagrams generating the quark-lepton 
mass matrices [||. Furthermore it is assumed that all the couplings, allowed by 
the AGUT gauge symmetries, of these Planck scale particles to each other and 
to the lighter Standard Model particles are of order unity. 
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In practice we construct the mass matrix elements between the known Stan- 
dard Model left- and right-handed Weyl fields one by one, by finding out first 
which quantum numbers need to be broken to make a given matrix element dif- 
ferent from zero. Knowing the needed breaking quantum numbers, we then ask 
which combination of Higgs fields in our model can together cause that breaking. 
Measuring everything in "fundamental" units we can say that, since everything 
is of order unity except the small VEVs, we just take the product of the VEVs 
for the needed combination of Higgs fields. Finally the mass matrix element 
in question is simply order of magnitudewise equal to this product. However 
we should mention that there is an ambiguity concerning which combination of 
Higgs fields to use, since the Abelian quantum numbers of our proposed Higgs 
fields W, T and ^ obey a Hnear relation: 3Qw — 9Qt — 2Q^ = (mod S quan- 
tum numbers). The resolution of this ambiguity is to choose the combination of 
Higgs field VEVs giving the largest value for the product. There is no problem 
in practice due to the large powers of T involved. 



3 Combinatorial corrections to mass matrices 

We will now discuss a modification of the fitting procedure we have used previ- 
ously j^]. The elements of the mass matrices are determined Q up to order of 
one factors to be a product of several Higgs vacuum expectation values (mea- 
sured in units of some fundamental scale Mp which we will take to be the Planck 
scale.) This is because the entries in the mass matrices come from Feynman 
diagrams such as that shown in Fig. [|. We have used such mass matrices to 
fit the fermion masses and mixing angles. In fact we have managed to produce 
fits with a smaller than would be expected, given the lack of knowledge of 
order of 1 factors. So it is sensible to try to make a more refined fit, which 
takes into account systematic factors, even though we will still have order of 1 
uncertainties which we can only average over. 

One effect we have ignored up till now are the combinatorial factors due to 
the internal ordering in the Feynman diagrams. That is, we have only taken 
a single Feynman diagram to calculate the order of magnitude of each mass 
matrix entry. What we should really do is sum up all Feynman diagrams with 
given initial and final states. If we restrict ourselves to tree level diagrams this 
means that we should consider the sum over all distinct orderings of the inter- 
actions with the Higgs VEVs. If we consider a diagram with order of magnitude 
(M^)°(T)''(^)°(S')'^($VKs), then the corresponding matrix element should now be 
calculated by summing over all '''^^''a^h\ci.d\^^''' diagrams. Since in our approach 
only the order of magnitude of each diagram is determined, we should add all 
these diagrams together with random phases. Therefore the order of magnitude 

of this matrix element is given by ^"+1^+"+^,+^^' ' (W)''{T)''{0'={S)'^{<i>ws)- 
It would seem at first sight that these combinatorial factors would greatly 



6 



¥ ¥ ¥ 



$w,9 \ W \ T 



bL Mf ^2 Mf ^3 bR 

Figure 1: Possible Feynman diagram for bottom quark mass. The crosses indi- 
cate the couphngs of the Higgs fields to the vacuum. Note that this is only 1 of 
6 possible tree level diagrams since the order of the interactions is arbitrary. 



modify our mass matrices and so vastly alter the fits we have already analysed 
without such factors. However, this is not quite true since the factors can to 
a large extent be absorbed into the Higgs VEVs. This means that we can get 
similar fits, with or without such factors, but with rescaled Higgs VEVs. But of 
course there will still be some differences due to these factors, since they cannot 
be exactly absorbed by redefining the Higgs VEVs. One important point is 
that now even if we assume that the Higgs field S has VEV (S) = 1, this 
field still affects the fit through its contribution to the combinatorial factors. 
So, whereas previously the quantum numbers of such a Higgs field could be 
freely absorbed by other Higgs fields without changing anything, we will now 
get slightly different fits depending on the precise definition of the quantum 
numbers of the Higgs fields. 

4 Mass matrices and Higgs quantum numbers 

The quantum numbers of the Higgs fields S, W, T, ^ and ^ws were origi- 
nally constructed by requiring them to give phenomenologically acceptable 
relations between the lepton and quark masses and their mixing angles. As 
we have already emphasized, the quantum numbers of the W , T, ^ and ^ws 
fields were only determined modulo those of the S field. We here arbitrarily 
select a standard set of such Higgs field quantum numbers, with "small" values, 
and parameterise the equivalent sets of quantum numbers by integer parame- 
ters (a, P, 7, S), defining how many S quantum numbers have been added. This 
leads to the Higgs fields having the following charges: 

Qs - -io,-l) (5) 
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Qw 


= ( g, 3,3, ^)+aQs 


(6) 


Qt 




(7) 


Q« 


= (^.-^-0,0)+70s 


(8) 


Qws 




(9) 



The orders of magnitude of the up quark, down quark and charged lepton 
mass matrices are given in terms of the Higgs fields as follows: 



Mu = Hu-P^s (10) 
Md = Hd^ws (11) 
Me = He^ws (12) 

where the order of magnitude of the Yukawa coupling matrix Hu is given by: 



(^l+c<-2/3+27+5i[^t2^2(-^t)2 ^2+a-2/3-7+5|^t j.2^ ^2+a-/3-7+5(-;^t-)2 j.^ > 
^l+c<-2/3+37+<5p^tr2^^t-)3 ^2+a-2/3+5pj/t j>2 g2a-l3+S (jy^)'^T 

the order of magnitude of the Yukawa coupling matrix Hd is given by: 



g-l-a+20-2j-5y^(j^r-j'2^2 ^-2-a+2/3-7-5p^(-2.t)2^ 5'2-3/3-7-<5y3^ \ 
5'-l-a+2/9-7-5yi7(-yt)2^ g-2-a+20-5-^Y (^^^.^^2 ' g2-30-Srp3 j 

^-2-2a+4/3-7-5|y2 ^2^1)4^ g-3-2a+4^3-S^r2^rpiy gl-a-|3-S^rJ' J 

and the order of magnitude of the Yukawa coupling matrix He is given by: 



^_l_a+2/3-27-5^(-yt)2^2 ^-2-a+2/3+37-<5^(-yt)2(-^t)3 g2-a-4^3+J-S^^J^4^^ 
g-l-a+20-5^-6^r^J^^(^2^5 c<-2-a+2l3-Sy^r(rpii-^2 g2-a-il)-2'y-5y^j.4^2 

5'-2-c<+5/3-37-5^^yt)5^3 ^l+2a-4/3-5^|yt^2y4 gi-a-pSy^j' 

The entries in the mass matrices represent the Higgs fields involved in the 
Fcynman diagram describing the trcc-lcvcl interaction relevant for each element. 
It is to be understood that the actual values of the mass matrix elements are 
given by the products of the expectation values of the Higgs fields involved. 
So, for example, (VK''')^ would mean {W)'^ and S*" means (S*)!"!. Also, for 
simplicity, we haven't included the combinatorial factors in the mass matrices. 
These factors must be included in the fit and are easily calculated from the 
powers of each Higgs VEV, as in the previous section. 
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So the charges of the Higgs fields used in our previous paper 0| were given by 
the choice a — 1, (3 — 1,7 = and i5 = — 1. However, without the combinatorial 
factors and with {S) = 1, the choice of a, /3, 7 and 6 does not affect the fit. But 
now we will include these factors and vary the integers a, f3, 7 and S to see what 
the effect on the fit is. Clearly large values of these variables will introduce large 
powers of (S) and large combinatorial factors which will change the character of 
the mass matrices. Since we have derived the model under the assumption that 
(S) ~ 1 and that the combinatorial factors are not very important, we should 
restrict ourselves to small values of the integer parameters a, j3, 7 and 5. Also 
this will satisfy our general requirement that the model should have "small" 
quantum numbers. We allow {S) to vary, in order to somewhat compensate 
for the combinatorial factors. We decided to allow values -1, and 1 for each 
parameter a, /?, 7 and 5. This gives us a total of 81 choices. We want to find 
the best fit among these 81, minimising a pseudo-chisquared function defined in 
Eq. ( [Tsl ) below. To do this we first made an approximate fit for all 81 choices, 
giving an average pseudo-chisquared x^^g = 2.7, and then chose those fits which 
had < 2.0. This left us with 14 fits, which we then analysed with higher 
accuracy. 



5 Results 

As in previous papers, we fit to the experimental values given in Table and 
use a pseudo-chisquared function to measure how good a fit is. Since we are 
making an order of magnitude fit, we use the definition: 

^ ^ .^^ ( "^exp ) 

where m are the fitted masses and mixing angles and mexp are the corresponding 
experimental values in Table ^. The was minimised by varying the Higgs 
VEVs, where in this paper we also vary {S) rather than fixing (5) — 1. The 
Yukawa matrices are calculated at the fundamental scale, which we take to be 
the Planck scale. We use the first order renormalisation group equations (RGEs) 
for the Standard Model to calculate the matrices at lower scales. Running 
masses are calculated in terms of the Yukawa couplings at 1 GeV. 

In this section we will comment on the fits in Tables |] and ^ in order to 
highlight the general features of different choices for the charges of the Higgs 
fields. We will start by making a comparison between using random complex 
0(1) factors in the fitting procedure or just random phase factors. 

As we have repeatedly stressed, we are only assuming some knowledge of 
the order of magnitude of the mass matrix elements. This means that, when 
calculating the eigenvalues of these matrices, we should at least consider each 
element to have an arbitrary phase. We do this by averaging over the calcu- 
lated eigenvalues, using many random choices of phases for all the elements. 



(13) 



9 



Table 1: Experimental values of masses and mixing angles used in the fits. All 
masses are running masses at 1 GeV except the top quark mass which is the 
pole mass. 



niu 4 MeV 

md 9 MeV 

me 0.5 McV 

Mt 180 GeV 

nib 6.3 GeV 

rur 1.78 GeV 



rric 1.4 GeV 
rris 200 MeV 
105 MeV 
K« 0.22 
Vcb 0.041 
Vub 0.0035 



Table 2: Results for fits without any 0(1) factors used. 



a 




7 


s 


(W) 


(T) 


(5) 


iO 


x' 


-1 


-1 


-1 


1 


0.0672 


0.0667 


0.487 


0.0331 


1.81 


-1 





-1 


1 


0.0857 


0.0522 


0.33 


0.0365 


1.34 


-1 





1 


-1 


0.0705 


0.0549 


0.720 


0.0422 


1.88 


-1 


1 


-1 


1 


0.0735 


0.0525 


0.686 


0.0331 


1.90 


-1 


1 


1 


1 


0.0857 


0.0498 


0.720 


0.0402 


1.68 








-1 





0.0610 


0.0851 


0.286 


0.0315 


1.47 








-1 


1 


0.0671 


0.0576 


0.464 


0.0402 


1.42 














0.0705 


0.0810 


0.259 


0.0365 


2.04 





1 


1 





0.0945 


0.0522 


0.876 


0.0365 


1.78 




-1 


-1 


1 


0.0740 


0.0576 


0.622 


0.0315 


1.91 




-1 


1 


1 


0.0777 


0.0548 


0.622 


0.0402 


2.03 







-1 





0.0777 


0.0575 


0.537 


0.0443 


1.87 







-1 


1 


0.0740 


0.0605 


0.512 


0.0331 


1.94 




1 


-1 





0.0992 


0.0548 


0.271 


0.0489 


1.68 
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Table 3: Results for fits including 0(1) factors. 



a 


a 


7 


5 


m 


(T) 


{S) 


(0 


X 


2 


-1 


-1 


-1 


1 


0.0741 


0.0635 


0.487 


0.0331 


1 


.57 


-1 





-1 


1 


0.0945 


0.0522 


0.347 


0.0331 


1 


.41 


-1 





1 


-1 


0.0857 


0.0522 


0.686 


0.0365 


1 


.59 


-1 


1 


-1 


1 


0.0894 


0.0525 


0.756 


0.0247 


1 


.26 


-1 


1 


1 


1 


0.0945 


0.0474 


0.653 


0.0365 


1 


.46 








-1 





0.0741 


0.0810 


0.286 


0.0300 


1 


.27 








-1 


1 


0.0857 


0.0548 


0.442 


0.0347 


1 


.40 














0.0816 


0.0735 


0.299 


0.0331 


1 


.37 





1 


1 





0.0945 


0.0522 


0.721 


0.0331 


1 


.62 


1 


-1 


-1 


1 


0.0857 


0.0522 


0.622 


0.0300 


1 


.44 


1 


-1 


1 


1 


0.0900 


0.0497 


0.622 


0.0383 


1 


.70 


1 





-1 





0.0900 


0.0522 


0.537 


0.0422 


1 


.79 


1 





-1 


1 


0.0816 


0.0549 


0.538 


0.0331 


1 


.64 


1 


1 


-1 





0.1042 


0.0522 


0.346 


0.0444 


1 


.85 



Table 4: Typical fit without averaging over 0(1) factors with a = — 1, /3 = 1, 
7 = 1 and 5 = 1. All masses are running masses at 1 GeV except the top quark 
mass which is the pole mass. 





Fitted 


Experimental 




4.2 MeV 


4 MeV 


nid 


4.7 MeV 


9 MeV 




0.98 MeV 


0.5 MeV 


nic 


1.22 GcV 


1.4 GcV 


rus 


340 MeV 


200 MeV 




83 MeV 


105 MeV 


Mt 


220 GeV 


180 GcV 


mb 


7.2 GeV 


6.3 GeV 


rrir 


1.17 GeV 


1.78 GeV 


Vus 


0.15 


0.22 


Vcb 


0.031 


0.041 


Vub 


0.0040 


0.0035 


Jcp 


9.4 X lO"*^ 


2 - 3.5 X 10"^ 




1.68 
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Table 5: Typical fit including averaging over 0(1) factors with a — —1, /? = 1, 
7 = 1 and i5 = 1. All masses are running masses at 1 GeV except the top quark 
mass which is the pole mass. 





Fitted 


Experimental 




3.1 MeV 


4 MeV 


TOd 


6.6 MeV 


9 MeV 


TOe 


0.76 MeV 


0.5 MeV 


TOc 


1.29 GeV 


1.4 GeV 


TOs 


390 MeV 


200 MeV 


TO^ 


85 MeV 


105 MeV 


Mt 


179 GeV 


180 GeV 


TUb 


7.8 GeV 


6.3 GeV 


m-r 


1.29 GeV 


1.78 GeV 


Vus 


0.21 


0.22 


Vob 


0.023 


0.041 


Vub 


0.0050 


0.0035 


Jcp 


1.04 X 10"^ 


2-3.5 X 10-5 




1.46 





This avoids any accidental cancellations between quantities of the same order 
of magnitude. However, we can also consider introducing random factors which 
are of order 1, since this will not change the order of magnitude of each element. 
Now, of course, we must decide more precisely what we mean by an 0(1) factor 
and how to average over the random variations. There is no unique way to 
decide how large a number should be before we no longer consider it to be of 
order 1. However, a reasonable choice is to say that a real number is of order 
1 if its natural logarithm lies between -1 and 1. So we will pick random 0(1) 
factors by taking the exponential of a number picked from a Gaussian distribu- 
tion of width 1. We will now compare the results from fitting after averaging 
over these 0(1) factors to the results without using any 0(1) factors (but still 
averaging over random phases.) Tables ^ and ^ give the values for the best 
14 fits with and without 0(1) factors. However, in order to illustrate in detail 
the differences, we have shown the fitted masses and mixing angles for a typical 
fit in Tables | and |. 

The most obvious differences between using 0(1) factors or not is in the 
ratio between up and down quark masses. The difference between these masses 
is always increased when 0(1) factors are used. The reason for this is that 
the down mass is produced by two different combinations of matrix elements. 
Since the down mass will typically be given by the root mean square average of 
these two combinations, random factors will generally increase this mass more 
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than other masses which are essentiahy only determined by one combination of 
matrix elements. In most of the fits this helps but, in a few, the ratio is already 
large enough without the 0(1) factors and this makes the fit worse. The 0(1) 
factors also seem to allow a better fit for the electron mass, which is too high 
without the 0(1) factors. This is because the down mass can still be fitted with 
a lower electron mass (and so lower up mass), if the ratio of down:up masses is 
increased. It can be argued that we could adjust the spread of 0(1) factors to 
tune this ratio. Then this would introduce an extra parameter in our fit, since 
there is no good way to fix it. However, we would argue that we have chosen a 
'natural' spread of these factors and that this spread can be considered, loosely 
speaking, as a definition of what we mean by 0(1). 

The other effects of using 0(1) factors are less important. In the third 
generation the 0(1) factors reduce the top mass and increase the bottom mass, 
while leaving the tau mass almost unchanged. However, these are small effects 
although this does usually worsen the predictions. In the second generation the 
0(1) factors increase the charm mass, which is still too low, but also increase 
the strange mass, which is always too high in our model. The muon mass 
is not changed much. These are again small effects and the improvement of 
the charm mass approximately cancels the worsening of the strange mass in 
the fit. However, this obviously increases the dominance of the contribution of 
the strange quark mass to the for the fits. The 0(1) factors increase Vus 
which is still predicted too small. They decrease Vet and this is usually bad but 
sometimes it is too high without 0(1) factors. Vub is increased, which sometimes 
helps, but often it is predicted too small without 0(1) factors and too large with 
0(1) factors so that, by chance, the 0(1) factors make almost no difference to 
the Vub contribution to the value of x^- Overall the 0(1) factors lead to a 
better prediction for the mixing angles and a small increase in CP violation. In 
general the 0(1) factors improve the fit, mainly due to the better fitting of the 
first generation masses and also the mixing angles Vus and Vub- Of the 14 fits, 
only 2 are worse with the 0(1) factors. In both cases this is because the down:up 
mass ratio is already large enough without 0(1) factors. For a = —1, (3 = 0, 
7 = — 1 and 6 = 1 the fit is only slightly worse with 0(1) factors and is still a 
good fit, since the worse up mass is compensated for by a better electron mass 
and Vus- For a = l,/3=l,7 = — 1 and 6 = many of the general comments 
do not apply. Here the down:up mass ratio is surprisingly smaller with 0(1) 
factors and the mixing angles Vus and Vub are almost unchanged. This leads 
to a worse fit with 0(1) factors, though we don't know why the down:up mass 
ratio changes in this way for this case. 

It is much harder to give any general comparison between the different 
choices of a, /3, 7 and 6- The values shown in Tables || and || have been chosen 
since we expected reasonably good fits from them. They are, in fact, the best 
fits out of all 81 combinations of a, (3, 7 and 6 with values -1, or 1. On general 
grounds we would expect worse fits from larger values of a, (3, 7 and 5. This is 
because large powers of the Higgs field {S) would be required for many matrix 
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elements. This would then suppress the elements unless (S) = 1. We would 
also expect large contributions from the factorial factors, but it would be un- 
likely that these could be balanced by the suppression due to (S) for all matrix 
elements. Therefore we would end up with some elements larger than expected 
and others smaller than expected. Since these effects have not been taken into 
account in our derivation of the model, they would almost certainly spoil our 
assumed relations. We will now just consider the cases with small values of a, 
P, 7 and 5 that have been displayed in the Tables ^ and ||. We have displayed 
the fitted for the fits. Tables ^ and || show the fitted masses and mixing 
angles for typical cases, with and without averaging over 0(1) factors. Also 
shown is the predicted Jarlskog CP violation invariant Jcp 

We can see that generally the largest contributions to come from the 
strange mass, up mass, electron mass, top mass and Vus- The strange mass is 
always fitted too high and is fitted the worst. It doesn't vary much among the 
fits, indicating that the fit minimises the strange mass since it is the largest 
single contribution to x^- Therefore the results clearly show that the strange 
mass cannot be accurately fitted within our model, no matter what values of a, 
P, 7 and S are chosen. So we are left with the conclusion that the strange mass 
must be accidentally small, if our model is to be believed. This is not necessarily 
a problem, since we only claim to fit order of magnitudes and a factor of 2 is 
not so large. 

We can see some trends in the way other masses and mixing angles are fitted 
with different values of a, /?, 7 and 6. There is generally a compromise between 
the up and electron masses (up too low and electron too high). The down mass 
is then largely determined by the ratio to the up mass, which varies between 
fits. This ratio depends mainly on 7 — certainly it seems that 7 = — 1 gives a 
better ratio than 7 = 1. This is presumably the main reason why most of the 
best fits (these 14 were the best out of 81) have 7 = — 1 and in fact the two very 
best fits both have 7 = — 1 . Unfortunately there doesn't seem to be any simple 
dependence on a, 7 or (5 for the predicted top mass or Vus- Since many of the 
fits are almost equally good, there is probably some 'unpredictable' dependence 
on the factorials which makes it harder to determine such a dependence. Since 
so many fits are fairly good, it seems that we have little chance of really choosing 
one scheme as the best. This is perhaps to be expected, since the choices of 
a, f3, 7 and S are in some sense just a fine-tuning of the general model. The 
model was 'derived' to include important relations between observed masses and 
mixing angles. Therefore, provided these variations don't spoil such assumed 
relations, we would not expect great differences in the overall fit. This is perhaps 
most obvious in that the strange mass is fitted 'equally badly' in the different 
variations of the model. 
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6 Making order unity concept more precise 



We shall now take into account how to define "order of unity" more precisely, 
when the number of particle species, or the dimension of some representation 
to which they belong, is so large that we do not want to consider it of "order 
unity". We shall, however, formally proceed in the estimate below as if the 
number of lines etc. in a Feynman diagram can be considered of "order unity" ; 
in particular we shall not consider the number of external lines to be more than 
of order unity. 

In the crudest approximation the following assumption should make sense: 

Assumption A: At the Planck scale all masses and couplings are 
of order unity and particles with all required quantum numbers — 
namely all quantum numbers possible — do exist. 

However we should like to work to a level of accuracy where, for instance, the 
number of colours Nc — which is say 3 — can no longer be counted as of order 
unity, but rather should be considered much bigger than unity. It is then not a 
priori obvious whether the above crude "order unity" approximation can even 
be made consistent to such an accuracy nor, more precisely, how it should be 
done. 

We here propose the following criterion for a successful more precise speci- 
fication of the order unity of couplings concept: 

The criterion: The main guideline we shall use is the requirement 
that composing two order unity amplitudes, say one N-point func- 
tion and one P-point function to one (N+P-2)-point function by 
combining via a propagator, should result in an order unity com- 
posite amplitude. This propagator is here assumed to only include 
particles that are not mass protected, so that they have Planck order 
of magnitude masses. 

The requirement that we only include propagators for particles which are not 
mass protected reflects the fact that, in assumption A, we "assumed the order 
unity to be valid for the particles at the Planck scale" . It means that the quarks 
or leptons themselves are not supposed to be included as propagator particles 
in this criterion. What we calculate are amplitudes due only to exchanges of 
heavy (Planck scale) particles, and all the quark or lepton propagators must 
be inserted explicitly later (after the use of the criterion and the construction 
of the Planck scale particle exchange amplitudes). An example of such a later 



insertion can be found in subsection 6.4 



Let us argue now for why we should take this criterion to be obeyed as a 
test for a successful implementation of assumption A: 

Consider an (N + P — 2)-point function. In tree diagram approximation 
we can ask for the contributions to this {N + P — 2)-point fimction coming 



15 



from diagrams with a propagator in a certain channel specified by, say, — 1 
of the external lines. For each combination oi N — 1 lines out of the N + 
P — 2 external lines wc have a contribution. The statement "particles with 
all required quantum numbers do exist" may be taken to suggest that none of 
these contributions are allowed to be much smaller than the other ones (i.e. we 
take them all to be of the same order of magnitude). Now such a contribution 
associated with a certain channel, in the sense of a selection of iV — 1 of the 
external lines, can be written as a P-point function connected by a propagator 
to an iV-point function. So such a single contribution is given by TV-point 
and P-point functions with a propagator connecting them. Since we take all 
the contributions equal in size, order of magnitudewise, and since there are in 
principle "only a few" (because wc assumed that the N etc. of the TV-point 
function is not large) they must all — order of magnitudewise — coincide with the 
full amplitude, the (A?'+P— 2)-point amplitude. But that means that we argued 
for the criterion to be fulfilled, in order that the assumption A can be said to 
be true. 

In the above argument we used the approximation — which is perhaps not 

so good in practice — that the number of external particles in the considered 
amplitudes were of order unity with sufficient accuracy. For the application 
to the quark and lepton mass matrices, this number is equal to the number 
of Higgs fields used plus 2 for the quarks or Icptons. However, we could say 
that the "factorial corrections" treated at length in section 3 — the main new 
feature in the fits presented in this article — were precisely to take into account 
factors arising from the number N for these amplitudes not being of order unity 
accurately enough. Actually it turned out that the effects we studied were the 



viVT factors, coming from including diagrams with the internal channels used 

being changed corresponding to the permutation of the external lines (or rather 
some of them). 

Wc shall now explain how reasonable "order unity" A/'-point amplitudes are 
proposed, and subsequently confronted with our criterion. First wc shall, for 
simplicity, discuss the case of there being no conservation laws, so that there 
are no forbidden couplings. 

6.1 Case of no conservation laws 

In order to define what the "order unity" should mean for the iV-point ampli- 
tude, it is convenient to consider Euclidean momenta, since then the poles in the 
amplitude arc avoided. The natural c;lioic;c; of definition is then that the TV-point 
function for an arbitrary set of particles, with Euclidean momenta of order unity 
in fundamental units, shall be of order unity in the simple mathematical sense: 



This is indeed what we shall propose, but there are a few things to have in 
mind here: 




A{pi,p2, ■■■,pn) ~ 1 (in Planck units). 



(14) 
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In usual notations, it is customary to use propagators for which the residues 
of "the" pole are normalized. But the propagator is the inverse — as a matrix — 
of the two-point function, and if we thus use our specification, Eq. (p^), to 
determine even the two-point functions, we loose the freedom to normalize the 
residues. 

Also we shall have in mind that there are many species of particles of the 
Planck scale mass type (we here do not include the mass protected quarks and 
leptons) and that the two-point functions make up a matrix, the inverse of which 



is the propagator, which also is a matrix. Our postulate, Eq. (14), means that 
the individual matrix elements in the two-point function are just of order unity. 
Thus the species, in terms of which we express the order unity amplitudes are 
NOT mass eigenstatesQ. We shall make these considerations of the order unity 
calculations more precise statistically, by assuming that we have a statistical 
ensemble of order unity coefhcients in the Lagrangian written in terms of fields 
for the various species, with a spread of order unity (in Planck units) so that 
they are practically always of order unity. Since the free terms — the kinetic 
term as well as the mass term — are then random, these species end up being 
rather random compared to the mass eigenstates of course. 

The propagator is then easily estimated to have matrix elements of the 
order l/(# of species). This may be seen by arguing that the eigenvalues of 
the inverse propagator are of the order -^/^ of species, while the eigenvalues 
of the propagator must of course then be of order 1 / yj (# of species). So the 
propagator is a factor of (# of species) smaller than the inverse propagator. 

Provided we use the propagator as the inverse of the two-point function 
matrix, it does not really matter if we multiplied all the 'N point functions by 
the A^th power of any number (e.g. the square root of the number of species 
•\/# of species). This is because such a factor for each line would be can- 
celled by the propagators, or in the case of external lines by the wave function 
(re)normalization. So actually we could equivalently have postulated the "order 
unity" rule as: 

A(pi,p2, ■••iPjv) ~ l/(# of species)^^/^^ (in Planck units) (15) 

This really may also be described as saying that the sum of the numerical squares 
of A over all combinations of external particles should be of the order of unity. 

But, before we accept one of these (essentially equivalent) ansatze, we should 
check that the criterion is fuUfilled: 

In order to check the criterion, we write down the double sumQ over the 
species at the end-points of the propagator connecting the A^-point and P- 
point functions to form a composite (N-|-P-2)-function. Now, since there are 

•^If wanted, the mass eigenstates would have to be found by dividing the mass term, as a 
matrix, by the kinetic term, also a matrix, and then diagonalizing the resulting matrix. 

*It should be remembered here that our species were neither mass nor kinetic term coeffi- 
cient eigenstates so that a DOUBLE sum is needed, the propagator being a matrix. 
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of species) terms in each of the two summations and we take each term to 
be a complex random number, these sums effectively each function as a factor 
of of species. Since the matrix elements of the propagator are of the order 
(l/# of species) in the notation of Eq. ( p^ ) and unity in the notation of Eq. (|l5|), 
we end up in both cases with a composite amplitude of the correct "order of 
unity" . 



6.2 Case with e.g. SU{Nc) symmetry 

In the case where there is a symmetry, hke colour symmetry under the group 
SU (Nc) say, we cannot simply use the rule for the case of no conservation laws 
given above, because a random N-point function would of course almost certainly 
violate the symmetry. Instead let us make the "order unity" assumption for N- 
point functions, by defining a probability distribution for amplitudes. Then we 
can say, in a precise sense, that an amplitude is of "order unity " when it is 
of a type that can be gotten with high probablity from this distribution. The 
probability distribution is then specified as we describe below. 

We first write down the cartesian product of the states of the N external 
particles counted as incoming (outgoing particles replaced by incoming antipar- 
ticles): 

\spl,cl> \sp2,c2> ■■■\spN,cN > e Hn- (16) 

Here the symbols "spl" ... "spiV" stand for the species (= "flavour" ) of the exter- 
nal particles, while "cl"..."cAf" denote their colour states. Often one imagines 
working with (external) particle states |spl,cl >, ...\spN,cN > that have, for 
instance, a definite weak isospin component along a certain quantization axis or 
definite values of some analogous colour Cartan algebra generators. However, 
for the purpose of formulating the statistical rules defining the "order of unity" , 
we shall here rather assume that the external states are random superpositions 
of various, e.g. colour and flavour, components. W.r.t. "flavours" (meant in the 
very general sense in which it is used here), it does not matter much to have 
this in mind, but for the conserved colour we risk to get many zeroes alternating 
with too big numbers for the amplitudes ^at, if we do not work with random 
states. So the following formulae are only meant to be valid for states that 
are random superpositions! An N-point function conserving colour is given by 
specifying a bra- vector < i?Ar|, belonging to the colour singlet subspace of the 



cartesian product Hilbert space Hn (see Eq. (16)), as the overlap 



AN{pi,spl;p2,sp2;...;pN,spN) = 
< Rn{pi,P2,-,Pn)\\spI,c1 > |sp2,c2 > • • • \spN,cN > . (17) 

We now make a random "order unity" N-point function by choosing this bra- 
vector < i?jv| randomly. We decide to take a "rotational invariant" direction 
distribution for it and a norm of the order of the square root of the dimension 
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of the space of singlets: 

II < i?jv|| ~ \/ dim{space of singlets) (18) 

This latter choice may a priori look rather arbitrary, but note that it is just 
the choice that corresponds to the overlap size of < Rn\ with any normalized 
singlet state \sing > in Hn becoming of order unity: 

< I < RN\sing > \ 

— averaging 

« 1 (19) 

Really this statistical definition of the order unity amplitudes is a bit super- 
fluous and we could more simply just state that the amplitude Apf becomes of 
the order of 



ANipi,spl;p2,sp2]...;pN,spN) = 
"0(1) - factor" < singletN\\spl, cl > \sp2, c2 > • • • \spN, cN > . (20) 



E 

singlet N 

where we sum over a basis set of singlet states. We should here comment on the 
normalisation of the states w.r.t. colour and flavour (i.e. species). We are using 
the normalisation of Eq. (|lj) for the species degree of freedom. However we 
have chosen to use the normalisation of Eq. (|l5|) for colour in the sense that, if 
colour conservation were broken without changing their order of magnitude, the 
amplitudes would be smeared out so that the colours would function as species 



with the normalisation of Eq. (15). In the case of no (colour) conservation 
laws discussed in subsection 6.1, the propagator with the notation of Eq. ( p!5| ) 
would be just of order unity. So it is expected that the colour part of the 
propagator should on the average be of order unity. Indeed in the case of a 
colour conservation law, as considered in this subsection, we can easily estimate 
the propagator as the inverse of the two-point function A2 . With external states 
belonging to a colour representation a of dimension n = dim{a^) and having 
colour components m and m' respectively, we obtain: 

TO, spN\A2\m', spP' >^ (21) 
~ (0{l)spN,spP' < 1 , 0\a, to; a, -m' >) (22) 
^ / 0(1) 

spN,spP 



dim(a) 



<5- - I (23) 



Propa-spN.spP'ipN)S^^, (25) 
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where then^ 

Jn 

Propa-spN,spP'{PN) ~ 0{l)spN.spP'-J]—^ — (26) 

Now we want to show that our order unity ampUtude definition, Eq. ( pO| ) , has 
the nice property of passing the test of being self-consistent under composition 
of an N-point function and a P-point function by propagator contraction. We 
write down the composite (N+P-2)-point function formed from the propagator 
contracted N-point and P-point functions as foUows: 

^ AN(pi,spl;p2,sp2;...;pN,spN) 

spN ,spP\a 

Propa-spN.,spP'(j>N)Ap{p[, spl'; ...■,p'p = pn, spP') (27) 

Here pi,p2, ■■■tPn~i,Pi,P2j ■■•,Pp_i are the four-momenta of the external par- 
ticles, and the first N-1 are attached to the N-point blob, while the last — and 
primed ones — are attached to the P-point blob. Furthermore sp stands for 
species or flavour with an analogous enumeration, and Propa-spN.spP'iPN) is 
the propagator, which is a matrix, Eq. (^), in flavour or species space. The 
momentum of the propagator, determined from four-momentum conservation, 
is pn — p'p- The triple sum runs over all pairs of species or flavours for each 
colour representation a in the sum. 

The "self-consistency" we want to show is that this expression, Eq. (p7|), 
equals AN+p-2{pi,spl\p2, sp2; ...■,pN-i, sp{N-l);p[, spV; ...■,p'p_^, sp{P -I)') 
order of magnitudewise, or rather that the statistical distributions are roughly 
equal for the two expressions. 

According to the property, Eq. (p^), the (N-|-P-2)-point function is of the 
form < Rn+p-2\3pI,cI > \sp2,c2 > ■ ■ ■ \sp{N - l),c{N - 1) > \spl',cl' > 
■ ■ ■ \sp{P — 1)', c{p — ly >, with < Rp+N~2\ having projection on average unity 
(order of magnitudewise) on all the singlet states \sing >. So we would be well 
on the way to proving the self-consistency condition, if we could show that the 
dimension of this singlet space were the same as that of the singlet space for the 
composed amplitude of Eq. ( |27| ) . So as an introduction to prove the remarkable 
"self-consistency" 

AN+p-2{pi,spl;p2, sp2; ...ipn-i, sp{N - l);p[, spV; ...;pp_i, sp{P - 1)') = 

An{pi, ■■■;pN,spN)Propa-spN,spP'{pN)Ap{p[, ...;pp ^pN.spP') (28) 

spN^spP' .a 

we remark that the number of colour contractions — i.e. of the singlets that can 
be used — is indeed the same for both sides of this equation. 

^Of course the order unity numbers 0{l)apjv,spP' depend on spN and spP' , but they are 
just of order unity. So we could logically just have left the indices out and written 0(1) for 
Oil)spN,spP' ■ 
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This can be seen in the foUowing way: Imagine the singlets "of the left hand 
side" constructed by first couphng the colour representations of the (A^ — 1) and 
of the (P— 1) external lines of the (7V + P — 2)-point function to respectively the 
a and a representations, finally collecting the singlets for the various a-choices. 
Classified this way it is easily seen that these singlets are in correspondence with 
"those of the right hand side" , since the a-classification corresponds to counting 
according to the representation a of the propagators. Thus we have the equation 
for the number of singlets: 

^singletsN+p-2 = i^singletspf{a) * ^singletsp{'a) (29) 

a 

where the representations a in the expressions ^singlets n{q^ and ^singlets p{a) 
just denote the representations to which the species, spN and spP' respectively, 
should belong. 

After having noticed the natural correspondence between the singlet com- 
ponents on the two sides of the to be proven equation, Eq. (|^, we go over to 
considering a single contribution — one of the ^singletsN+p-2 in Eq- ( p9|) — to 
An+p-2- a contribution from a single "jsmgZei >" out of the ^singletsN+p-2, 
called say s, to Ajv+p_2 (the s-given representation a in the channel between 
the P-1 and the N-1 has dim{a) — n ) is 

^N+P-2 — 

< singlet{s)\spl,pi, cl > \sp2,p2, c2 > • • • \sp{N — l),pjv-i, c(7V 
•|spl',p'i,cl' > \sp2',p'^,c2' > ■■■ \sp{P-iy,p'p_,,c{P-iy > 

= X! A%\spN,pN,m > Propa-spN.spP'iPN) 
rn,spN ,spP' 

Ap(with p'p outgoing) I spP', p'p, TO > 
= ^ A%\rn > Propa-,spN,spP' {PN)VnAp\ - m > 

TTL.SpN ,SpP' 

~ X! A%\rn > Propa-spN,spP'{PN)Ap\m' > 

rh,m' ,spN,spP' 

The |to > and |m' > symbols denote that the iVth, respectively Pth, one of 
the external particles on the part of amplitude symbol A^|m >, respectively 
Ap\m' >, is a particle with colour component to, respectively to', while the 
irreducible representation of the colour degree of freedom is given as a conse- 
quence of the s-index. The first step, Eq. (^l|), consists in artificially writing 
the contribution A|^_|_p_2 from a given singlet s as products of amplitudes 
and Ap, with a completeness sum over all colour states. We use an orthonormal 
basis |to >, namely \rn >< ml ~ 1. This is consistent with the order of 
unity flavour normalisation, Eq. (|l4|), of the amplitudes ^a'+p-2i ^p- 



-1)>. 
(30) 



(31) 
(32) 

(33) 
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As a further elaboration, we introduce an artificial double summation for the 
flavour (species) indices spN and spP' over # of species values. The statistical 
increase by a V# of species-factor for each of the two summations is compen- 
sated by inserting the propagator, Eq. (p6|). This is just repeating the result. 



obtained in subsection 6.1, that the number of species is irrelevant for splitting 
an amplitude Aj^^p_2 into amplitudes A^^ and Ap connected by a propagator. 

We remark that the amplitude denoted Ap(with p'p outgoing) | spP' , pp , rfi > 
in Eq. ( ^T] ) is not normalised according to the same convention as the "order 
unity" amplitude obtained by just crossing our amplitude Ap. This is because 
we do not normalise the state with outgoing pp according to our "order unity" 
rule. We rather use the relation 

< 1, 0|a, m;a, -m >w — £ii]=. (34) 
^y dim{a) 

to obtain: 

Ap(with p'p out going) I spP',pp, TO >« ^ dim{a)Ap\ — to > (35) 

The reason for this normalisation is that Ap(with p'p out going) |spP', p'p, to > 
then appears as a subpart of the Clebsch-Gordan construction of the Ajv+p_2 
amplitude in the simplest way, i.e. without any extra dim{a) factors. We are 
really just using Eq. (|3j) to convert the notation from outgoing to ingoing. So, 
in Eq. (^2|), it is meant that both the \rn > and | — to > states are counted as 
ingoing. To get to the result Eq. (^), we use the following rule — which we often 
use in these random order one treatments: we take it that the sum of a series of 
statistically roughly independent terms of the same size equals a typical term 
times the square root of the number of terms. 

From Eq. (^3|) we get by summation over the different contributions, i.e. by 
summation over s, that 

An+p-2= ^ Propa-spN,spP'{pN)AN\m > Ap\m' > . (36) 

But this is just the self-consistency equation we wanted to test! Therefore the 
order of magnitude unity, as we defined it for the (N-t-P-2)-point function, co- 
incides with what you get by composing an N-point function with a P-point 
function, that are each of order unity in our sense, with an intermediate prop- 
agator (that is of course summed over) . This fact we take as very suggestive 
for our chosen definition to be a good one. It means that it passed the 'self- 
consistency check', the criterion we gave at the beginning of section 6. 

6.3 Extraction of definite colouri (i=l,2,3) states 

For the applications of the order unity amplitudes, it is not so nice to only have 
them defined for "random" colour states. We think of states with respect to 
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the three different SU (3) groups in our AGUT-model, when we here talk about 
colour states. However, for simplicity, in this section we write formulae for only 
one of the "generation corresponding" colours. It would of course be useful to 
extract expressions for amplitudes of order one for states with definite colour 
indices. 

As an example we shall take the cases of amplitude contributions that can 
occur in the large limit, where we cannot use the e-symbol in colours because 
it has (infinitely) many indices and only a finite number of external particles. 
Really let us, for our example here, imagine that to the singlet-contributions 
correspond amplitudes that can be simply written by means of colour kronecker 
deltas in the fundamental representations (= triplets for SU(3)) Sap^ where then 
a and (3 run over the Nc colours. That is to say we take the contribution from 
one of the < singlet\ states to be proportional to a product of such kronecker 
deltas Sap 

< singlet\spl, cl > \sp2, c2 > • • • \spN, cN > 

= "normalization" ^J^<5q^^^^ • product of the cl"' etc. (37) 

where there can then be various numbers of indices attached to the various 
external states depending on their representations. Here the "normalization'^ 
factor is to be calculated so as to get the random states amplitude to come out 
right. If we only care for getting the large N^, factors right, we would not need 
to care for whether we symmetrise states with a few indices. Furthermore, in 
the cases of interest in our model, we have for one generation of colour all the 
time only triplets, singlets or antitriplets. So let us, as the example, consider 
that we only have contractions of triplets between the various external particles. 

With < singlet\singlet >= 1 the statistical average of the square modulus 
of the expression in Eq. ( ^7|) will be 

Averaged < singlet\spl, cl > \sp2, c2 > • • • \spN, cN > p) = 

nin2 • • • riAT 

(38) 

where is the dimension of the representation (of colour) to which the random 
state \spi,ci > belongs. Or in other words 

< singletlspl, cl > \sp2, c2 > • • • IspN, cN >w ^ — . (39) 

y/nin2 ■■■nw 

Now putting the random states (in colour space) into the right-hand side of 
Eq. (^ leads to the result 

"normalization" ^JJ^^Oift^ ■ product of the cl"' etc. 

« "normalization" —^^=^^= (40) 
^nin2 ■■ - nN 
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in the case we considered. We here used that two random triplets/iVc-plets, c2" 
and say, have components of the order of and the rule for summation of 
equal and independent terms giving a factor of the square root of their number. 
For each Saf3 there are only Nc non-zero terms to sum and thus only one factor 
i/jVc- Hence we have, for example 

So we conclude that, for the consistency of Eq. (^^, wc must take 

"normalization" « -—^^=^^= (42) 
^nin2 ■■■UN 

without any further Nc factors. 

Now the full "of order unity" amplitude was defined in Eq. (|2^) as the sum 
over all the possible singlets. In cases where we only have singlets that can be 
written as simple kronecker delta contractions, as in the "large Nc" case we 
treated here, it follows that the full amplitude is obtained by summing up — one 
for each singlet — the expressions 

< singlet\spl, cl > \sp2, c2 > • • • \spN, cN > 
"0(l)-factor" 



]^'5c(i/3i ) • product of the cl^'etc. (43) 



^nin2 ■■■UN 
to get 

An{sp1,c1]sp2,c2; ■ ■■■,spN, cN) 

y '' Oj^y^^^^O''' (USc fj]- product of the cl"»etc. . (44) 
the smglets \ t / 

The fourth root normalisation factor could be distributed through Eq. ( ^3| ) 
or (^) so as to give a factor following each kronecker delta. In other words 
if, instead of the kronecker deltas above, we use the combination 

^"'^ (45) 



/Nc 



where Nc is the number of colours (for the S'C/(3)-group of the generation in 
question, of course Nc = 3), we do not need the fourth root denominators. That 
is to say we can write e.g. Eq. ( ^3|) as 

< singlet\spl, cl > \sp2, c2 > ■ ■ ■ \spN, cN > 

= "0(l)-factor" • product of the cl"'etc. (46) 
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6.4 Implications of our order unity choice 

Now that we have chosen a prescription for making the order unity concept 
more precise, what does this choice imply for the possibihty of finding large 
factors in our mass predictions? 

As an example we consider the characteristic prediction of our model: that if 
the diagonal elements dominate the mass matrices, then the quarks and charged 
lepton in the same generation have order of magnitudewise the same masses. 
Could this prediction be changed by, say, an factor? 

In the model we use we have assumed that there are several Dirac fermions, 
as well as other particles, with masses of the order of the Planck scale for all 
quantum number combinations we can think of. So, in particular, the quantum 
numbers of say a right-handed b-quark occur on several particles which are 
very heavy. This should be understood in the terminology of Wcyl particles as 
follows: There are several right-handed Weyl fields with the quantum number 
combination of the right-handed b-quark and in addition several, but one less, 
left-handed Wcyl fields with this quantum number combination. So just one 
linear combination of these right-handed Weyl fermions is left without a partner 
and is thus massless compared to the Planck scale. It only gets its mass, at the 
end, by being paired with a Wcyl field having the quantum numbers of the 
left-handed b-quark and, only then, under the influence of the Weinberg-Salam 
Higgs field and other Higgs fields in our model. It thus gets a mass much smaller 
than the Planck scale and this particle is the b-quark. 

We imagine now that we have integrated out the fields of mass of the order 
of the Planck scale and replaced their eff'ects by effective N-point amplitudes for 
the lighter particles. For example the cfi^cctivc amplitude {l,r,W,T, ...) 

that describes the scattering amplitude between the various particles described 
by the symbols Z, r, W, T, ... which stand for their names and states: the I and 
the r arc the left- and right-handed quarks or Icptons. and the W, T, ... stand 
for our various Higgs fields. We shall be especially interested in the case of 
the states of the Higgs fields W, T, ...being those superpositions in which they 
occur in the vacuum. But we just wrote "W, T, ..." as an example and we can, 
of course, write a similar expression with any combination of our Higgs fields. 
Further h is the number of the Higgs fields W, T, etc. and the upper index 
''{massy^ means that wc ignore the kinetic part, or equivalently that we take 
the fermion external momenta to zero. 

If the external states W, T,... are thought of as the vacuum condensate 
states - having of course zero four momentum the A^^'^'^^'^ (l, r, W, T, ...)-amplitude 
becomes the two point function corresponding to the efi'ective mass term in the 
Lagrangian for the Dirac particle formed from the Weyl particles I and r. 

If the kinetic energy term in the Lagrangian formulation were simply ^ times 
a unit matrix, we could see that this "two" -point amplitude ^2+^**^ mass- 
matrix (in colour space), but now the kinetic term coefficients are also random 
corresponding to being of order unity in the sense which we defined. We could 
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renormalize the fields so as to bring the kinetic term to the standard form 
or, equivalently, we could divide by the ratio of the kinetic term to j^. Since 
the propagators propR and propL are defined as inverses of the kinetic terms, 
we really have to correct the mass-amplitude Aj™^**"^ by a factor . We now 
consider the mass squared for the Dirac particle formed, which can be gotten 
from the relation 



Tr 



4™r^(/,r,w^,r,...)^^4+r^(r,;,w^,r,...)^^ 



dim (representation) m^. (47) 



Here these propagators are really just matrices approximately proportional to 
the identity « V^'^m m' colour space, where n is the number of colour 

states for the repesentation of ? or r (assumed equal), n = cizTO(representation). 
We assume that I and r belong to the same representation with respect to colour 
in order that the mass term shall not violate colour conservation. The trace, 
Tr, stands for the trace in the space of various-really colour-states in the I and 
r channels. If, by colour symmetry, they all have the same mass (squared) we 
just get the dimension of the representation of the I or r state spaces from this 
trace-summation, explaining the (iirn(representation) factor on the right-hand 
side of Eq. (B^). 



From our prescription Eq. (20) for "order unity" amplitudes, we have 



A2+h{r;l;W;...;T) = 
''0{1)- factor'' <singletN\\r>\l>\W >■■ -IT > . (48) 

singlet jy 

provided the external states are taken as "random" states — random colour su- 
perpositions first of all. Now, however, we want to consider definite colour states 
and we can make use of the estimate, Eq. (^), in subsection 6.3 to obtain: 

A2+hir;l;W;..-T) (49) 

= "O(l)-factor" i Yl j • product of the cl"*etc. (50) 

the singlets ^ ' ^ 

We have to imagine that, for the various Higgs fields W, T, etc., we have some 
colour contraction kronecker deltas between different generation colours. We 
note that all the Higgs fields in our model must conserve the diagonal colour 
group, which is identified with the QCD-colour group. So the typical Higgs field, 
like T in our model, must take expectation values which are invariant under 
this diagonal SU{3), although they break spontaneously the separate SU{3)'s 
of the various generations. That is to say that, for example, T could be triplet 
and anti-triplet under the 2nd and 3rd generation S'C/(3)'s respectively, and 
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thus its vacuum expectation value would be proportional to an inter-generation 
kronecker delta, having one index of generation 2 and the other of generation 3. 

In the simple case where we can neglect contributions containing e-symbols — 
the large Nc approximation — we would simply get several kronecker deltas con- 
tracted in circular chains, because at the end there should be no free indices 
in a colourless expression such as Eq. (^7|). Each ring of circularly contracted 
kronecker deltas is quickly reduced to just the trace of one kronecker delta and 
thus to Nc, the number of colours, which is the same for all the generation colour 
groups. The dominant term will be the one corresponding to that \singlet > 
which leads to the highest number of circular chains, because that will give the 
highest number of A*"c-f actors. 

For example for a diagonal mass matrix element, in which we have the same 
generation of colour index on the I and the r states, it will pay best to get 
a kronecker delta directly contract / and r; otherwise potential further con- 
traction loops would miss the chance of existing, because the Higgs fields put 
between the I and the r in the contraction chains would already have been used. 
We can use this observation to argue that the trace on the left-hand side of 
Eq. ( p7| ) for a quark of a certain generation will be just bigger, by having one 
contraction-loop more, than the corresponding lepton in the same generation. 
That means then that this trace will be Nc times bigger for a quark than for a 
lepton in the same generation — if their masses are given by the diagonal mass 
matrix elements (i.e. except for c and t in our model). From the presence of 
the factor dim (representation) in Eq. (|4^), it then follows that to this "large 
Nc approximation" , when diagonal matrix element dominate, the quark and 
charged lepton masses in the same generation are equally big (up to order unity 
factors which do not contain Nc factors, as one could perhaps have feared.). 

We have thus seen that one of the major predictions of our model — the order 
of magnitude degeneracy of masses in the same generation as long as they are 
dominated by diagonal matrix elements — is NOT modified by any big factor 
even if the colour number, Nc = 3, should be considered "bigger than of order 
unity" . It follows that one can expect a higher accuracy than having to ignore 
a deviation by a factor A^c = 3 in our predictions. Indeed our fits are accurate 
to better than to a factor 3. We only showed this for the diagonal mass matrix 
elements and with ignored e's, but it is not difficult to extend the argument to 
include the possibility of e's for the intrageneration mass ratios. 

It is somewhat complicated to calculate the colour-counting corrections for 
other matrix elements, but preliminary studies suggest that even for the off- 
diagonal elements the corrections to our fit are not very big. 

7 The expected uncertainty 

According to our discussion of the many different permutations of the order of 
the Higgs field tadpoles in the Feynman diagrams contributing to a given mass 
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matrix element, it could easily turn out that the mass matrix elements appear in 
our model as sums over rather many terms, all having in principle the same order 
of magnitude. There is also the possibility of having several fundamental scale 
particles with the same quantum numbers. If indeed there are many terms of 
the same approximate size, adding up with some random order of unity factors, 
we should be able to apply the central limit theorem to conclude that the sum, 
considered as a stochastic variable z, will have a Gaussian distribution in the 
complex plane (really a two-dimensional Gaussian distribution). 

But now we mainly care for the order of magnitudes of the mass matrix ele- 
ments and calculate the x^-like quantity for the logarithm of the variable rather 
than for the variable itself. So we shall now consider what the spread in the 
logarithm is for a variable that is Gaussian distributed. We shall assume, what 
because of their assumed random phases is true for the mass matrix elements 
in our model, that the average is zero. Then, for dimensional reasons, there 
is no way in which this fluctuation can depend on the width of the Gaussian 
distribution of the quantity itself. 

We now calculate the root mean square fluctuation of the logarithm of the 
stochastic variable z: 

< (logz)2- < logz >2> (51) 



/ exp(-|z|")(log|z|)^d"z _ / / cxp(-|z|")log|z|d'^z 

/ OXp(-|z|2) ^ J OXp(-|2|2)d2^ 

J^°° cxp(-t)(l/2*logt)^(it / JJ^ exp(-t)l/2*log tdt 

cxp(-'f)(it \^ JJ^ exp(-t)dt 



2 



(52) 
(53) 

- = ^ = 0.4112 (54) 

If all the masses in nature really followed our model in a statistical way they 
would deviate in the logarithm from the predicted values statistically by a factor 
of •\/'0'(l)/4 = 0.6412. Now we make a fit with 12 - 4 = 8 degrees of freedom, 
and so we expect the pseudo-chisquared of Eq. ( [T^ ) to be 8 x V''(l)/4 = 8 x 
0.4112 = 3.29. This to be expected pseudo-chisquared is of a very similar order 
of magnitude to the pseudo-chisquareds which we find in our fits. So roughly 
our agreement is similar to the expected one. But actually we got a somewhat 
better agreement than one should have expected for our presented fits! 

We should bear in mind, however, that the fluctuation of the pseudo-chisquared 
itself is of the order of 35%, so that the predicted pseudo-chisquared of 3.3 ± 1.2 
should be counted as only one standard deviation off if it were measured to be 
2.1. It is, of course, clear that if we fit 81 models and look at the ones with the 
smallest pseudo-chisquared, we should get some that have an anomalously low 
pseudo-chisquared. However, we investigated the distribution of our original 81 
approximate fits and found a pseudo-chisquared fluctuating around 2.1 with a 
spread of 0.9. This means that the major part of these fits work too well! A 
possible explanation for this is that all these 81 fits are not really very different 
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and should be considered roughly as just one fit, accidentally fitting too well by 
one standard deviation. An alternative explanation could be that some of our 
predictions are actually valid more accurately than in our model due, say, to 
some physics which our model has missed. For example we have the order of 
magnitude-but not exact-SU(5) mass predictions in our model. It could easily 
be that an exact SU(5) could lead to too good fits for our model. 

8 Conclusion 

We have improved our previous AGUT model fits ^ to the quark-lepton mass 
spectrum in a couple of ways: 

First we have taken into account the number of different permutations of 
the Higgs field attachments to the Feynman diagram contributing to a given 
mass matrix element. This is a necessary correction that should be expected 
to be there. Secondly, in order to readjust the fit, it seemed necessary to let 
the VEV of the Higgs field S, which was previously set equal to unity, also be 
a fitted parameter. Because of the S field now being also relevant, there are 
many variants of the model obtained by different choices of the gauge quantum 
numbers of the Higgs fields W, T, etc., which were previously only determined 
modulo those of S. This gives a lot of models, but they are in reality not very 
different. 

We proposed a more careful definition of what an order one coupling assump- 
tion shall mean and investigated whether it could give rise to correction factors 
of order Nc = 3. Such large corrections do not in fact appear, at least in our pre- 
diction of intrageneration mass degeneracy (except for the top and charm quark 
masses, which are dominated by off-diagonal mass matrix elements). Thus there 
appears to be no reason why we should not trust our naive calculations to an 
accuracy better than a factor of 3. 

So, as far as the charged quark and lepton mass matrices (neutrinos are a 
problem for several reasons are concerned, we managed to successfully fit 
the 9 masses, 3 mixing angles, and the CP-violation strength with our AGUT- 
model, i.e. with the gauge group SMG^ x U{l)f, and five Higgs fields W, 
T, ^, S and (pws- In addition to the usual Weinberg Salam Higgs vacuum 
expectation value — which in our model is replaced by 4>ws — we only used the 
vacuum expectation values of the four other Higgs fields as free parameters in 
our fits. It should then be remarked though that: 

1) With the slightly complicated gauge group representations of our Higgs 
fields, some fitting of discrete quantum numbers could be said to have taken 
place. 

2) The expectation value of the Higgs field S is so close to unity that we 
only barely need it in the fits, and so we should hardly count its value as a 
parameter. Its value is seldom fit to be less than 0.3 and it never needs to be 
more than 1. 
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3) Really we have a significant amount of ambiguity in the discrete quantum 
numbers mentioned under 2) in the sense that many proposals, deviating from 
each other by adding the quantum numbers of the field S to one or several of 
the other Higgs fields, give very good fits. 

The fits are so good that many of the fitting proposals mentioned under point 
3) even fit the data slightly better than could be expected. We calculated an 
expectation for the , defined in Eq. ( |l3| ) for the logarithmic comparison with 
data, from the hypothesis that the mass matrix elements at first have Gaussian 
distributions. This expected value for is 3.3 ± 1.2, which roughly agrees with 
the values of our fits. For a lot of the 81 models considered, we actually 
obtained a somewhat smaller x^ than this theoretically predicted one. So we 
even have reasonably good agreement concerning the degree of accuracy of our 
fit and must conclude that the model fits the charged fermion mass matrices at 
least as well as one should expect. 
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